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Spin waves, or magnons, in the anti-ferromagnetic a phase of solid oxygen provide a novel mecha- 
nism for ultra-cold neutron (UCN) production. Magnons dominate the energy exchange mechanisms 
for cold neutrons and UCN in solid a-oxygen, much in the same way as do phonons in solid deu- 
terium superthermal UCN sources. We present calculations of UCN production and upscattering 
rates in S-O2. The results indicate that S-O2 is potentially a much more efficient UCN source 
material than solid deuterium. 
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I. INTRODUCTION 

Ultra-cold neutrons (UCN) provide a low energy sys- 
tem in which to carry out fundamental neutron physics 
experiments to a much higher degree of precision than 
that previously achieved with cold neutron experiments. 
UCN typically refer to neutrons with velocities less than 
8 m/s. Potential experiments using UCN include the 
measurements of neutron /3-decay parameters to test the 
Standard Modelji^ the search for the neutron electric 
dipole moment to test time reversal symmetry;^ and 
the search for neutron-anti-neutron oscillations to test 
baryon number conservation. To date, some of these ex- 
periments have already been conducted with UCN but 
most of them are still statistically limited. To increase 
the UCN flux, Golub and Pendlebury 6 proposed a "su- 
perthermal" source, in which the UCN production rate 
is enhanced through the efficient inelastic coupling of 
cold neutrons to the phonon modes of a condensed mat- 
ter sample. The UCN loss rate is suppressed in these 
sources by the reduction of the thermal phonon popu- 
lation. Superthermal UCN sources have been success- 
fully realized in super-fluid helium 7 and solid deuterium. 8 
Here, we propose solid oxygen (S-O2) as a superthermal 
source, which shows considerable promise for enhanced 
UCN production via a novel mechanism. 

Solid cv- 16 02 has manifest magnetic properties. Its 
ability to interact with cold neutrons and UCN, through 
neutron-magnon inelastic scattering, provides potential 
advantages which stem from the energy dispersion spec- 
trum of its anti- ferromagnetic magnons. The large 
neutron-magnon inelastic scattering rate together with 
oxygen's small nuclear absorption cross section suggest 
that a S-O2 UCN source may well compete favorably with 
solid deuterium UCN sources. 

This report is organized as follows. In section II, we 
calculate the neutron-magnon and neutron-phonon scat- 
tering cross sections in solid oxygen, after reviewing the 
relevant magnetic properties and magnon dispersion re- 
lations. In section III we compare the predicted perfor- 
mance of a solid oxygen UCN source with the known 
performance of typical solid deuterium sources, and dis- 
cuss possible mechanisms for UCN loss. We summarize 
our results in section IV. 



II. NEUTRON-MAGNON SCATTERING 



A. Magnetic Scattering Length 



In solid oxygen at low temperatures (T<23.9K), 
thermal fluctuations quench and long range anti- 
ferromagnetic ordering is established. This results from 
a slight distortion of the lattice structure from a rhom- 
bohcdral (/3 phase) to a monoclinic symmetry (a phase) . 
The resulting 2-D anti-ferromagnetic structure exhibits 
spin wave excitationsAiS For our work on UCN, we will 
be concerned with S-O2 in the a-phase below 10K, where 
superthermal UCN production is possible. 

The magnetism of solid oxygen originates from the un- 
paired 2p electrons in the oxygen atom, two of which align 
and form a molecular ground state E~ with a net spin 
5=1 and an orbital angular momentum I = 0. For 
a diatomic molecule made of two 16 atoms (naturally 
abundant), the null nuclear spin and antisymmetric elec- 
tronic wave-function lead to a selection rule allowing only 
rotational states with odd parity»ii The absence of half 
of the rotational states is an advantage in the application 
of S-O2 as a UCN source, as will be discussed later. 

A neutron interacts magnetically with an oxygen 
molecule. The interaction potential can be described as 



V = n ■ H = ii ■ V 
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where /1 = 7/iArcr is the neutron magnetic moment with 
neutron spin ex, and fie — ~^Hbs is the electron magnetic 
moment with electron spin s. This magnetic coupling, 
which is the origin of neutron-magnon scattering, has an 
intrinsic scattering length of 5.4 frnii^ This is comparable 
to typical neutron- nuclear scattering lengths, involved in 
neutron-phonon scatterings. 

The differential cross section for neutron scattering in 



the presence of the potential (JTJ can be written: 
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Here hi is the initial neutron momentum vector, oriented 
relative to a crystal axis by a solid angle fi*, flf is the 
solid angle for final neutron momentum kf relative to 
the same crystal axis, and Aj, A/ are the initial and final 
states of the oxygen molecule. The above formula aver- 
ages over the initial states of the target molecule and the 
neutron spin. From Eq.ljl). the interaction matrix ele- 
ment is reduced to a scalar product between the neutron 
spin and the coherent lattice spin form factor Q, which is 
defined using the neutron momentum transfer k = kf — hi 
and the electron spin §i as 



Q = y^e tK ' ri K x (sj x fi). 



(3) 




FIG. 1: A schematic of the lattice structure of solid Q-O2. 



This 2-D anti-ferromagnetic system can be described 
by the Heisenberg hamiltonian 



H = -2 
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■J2(-D S 2 xi -D's 2 yi + D's 2 zi ), (5) 



The differential cross-section is conventionally simplified 
as a product of the square of the scattering length, the 
kinematic factor -^- and the scattering law S(«, oj)m^ The 
effective scattering length ro, characterizing the strength 
of magnetic neutron scattering, is given by 
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where 
and 7= 



— = 2.82/to is the the classical electron radius 
-1.91 is the neutron gyro-magnetic ratio. The 



detailed physics of neutron interactions in solid O2 is ab- 
sorbed in the scattering law S(k,lu). 



B. Anti-ferromagnetic Magnons in S-O2 

Calculation of the scattering law requires knowledge of 
the energy dispersion relation in the magnetic lattice of 
S-O2. As shown in Fig[fl solid a-02 has a primitive unit 
cell with a distorted hexagon geometry^ There are 14 
oxygen molecules in the unit cell. On the basal plane, the 
b-axis is defined to be along the magnetization direction 
(which is also along one side of the hexagon), and the 
a-axis is assigned to be orthogonal to the b-axis. Due 
to a larger lattice constant along the c-axis, couplings 
between magnetization planes were found to be negligible 
compared to the in-plane couplings. The dominant 2- 
D magnetic structure are comprised of two interlacing 
sub-lattices, each carrying the opposite magnetization. 
Every lattice site is surrounded by 4 nearest neighbors 
on the opposite magnetic sub-lattice, and 2 next-nearest 
neighbors on the same magnetic sub-lattice. 



where Jy characterizes the strength of exchange interac- 
tions between molecular spin s$ on site i and Sj on site 
j. D represents the self-energy of the longitudinal spin, 
and D' the self-energy of the transverse spin (here, the 
b- and a-axes are parallel to the x and y-axes, respec- 
tively). Note that even though the out-of-plane coupling 
is weak, the self-energy term D's 2 alone in the hamilto- 
nian can result in significant contributions to scatterings 
in the z-direction. 

The hamiltonian of Eq. (JSJ) leads to the dispersion re- 
lation of magnon excitations 9 : 



u> a (q) = ^{2J ~2J^-2J^ 



D) 2 -[2JL+{-l) a D>] 2 . 

(6) 

Here u) a (q) = W g , a is the energy of a magnon, with mo- 
mentum q. The angular dependence of the dispersion re- 
lation comes from J$ and JL, which are the Fourier trans- 



form of the coupling coefficients J^ 
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fined on the sites on one sub-lattice, and J' defined in the 
same way but on the opposite sub-lattice. Following the 
analysis of Stephens and Majkrzak^ we take the numer- 
ical values of the spin coupling Jnn between the nearest 
neighbors to be -2.44 meV, the value Jnnn between the 
next nearest neighbors to be -1.22 meV, and J± for the 
out-of-plane (z direction) coupling to be zero. Also, we 
use L>=0.132 meV and D'=0.118 meV. Here, a = {0, 1} 
is the branch index of magnon excitations, which include 
one acoustic and one optical mode. 

Fig. illustrates the dispersion relation at several dif- 
ferent neutron scattering angles, together with the free 
neutron energy E n = h 2 K 2 /2m n . For certain values of 
the incident neutron momentum, incident neutrons can 



fa)E 




(b)E 



c E 




fd)E 



FIG. 2: The dispersion relation (Eq.©) of magnons in solid 
a-02 plotted for momentum transfer along (1,0,0) in (a), 
(0,1,0) in (b), (0,0,1) in (c), and (1,1,1) in (d). There are two 
branches (acoustic and optical). In (a,b,c) different periodici- 
ties originate from different lattice constants on these crystal 
axes. In (d), momentum transfers along a non-symmetry axis 
produce a non-periodic dispersion relation. Intersections with 
the parabolic free neutron dispersion set the four momentum 
(energy, momentum) for UCN production. 



lose essentially all their kinetic energy in the excitation 
of single magnon, and emerge as (or "downscatter" to) 
UCN. These values are determined by the points where 
the free neutron curve intersects the magnon dispersion 
curves in Fig. [5] Note that this process is coherent, in- 
volving the simultaneous conservations of energy and mo- 
mentum in the scattering process, and is analogous to 
the creation of UCN via phonon production in superfluid 
4 Hc. We are only interested in the scattering conditions 
in which UCN are produced in the final state. Processes 
involving the creation of multiple magnons are also pos- 
sible, but these are higher order effects which we expect 
to result in corrections on the level of a few percent and 
we neglect them in our calculation. 

The energies of typical excitations in condensed matter 
systems range from tens to hundreds of Kelvin. This is 
comparable to the kinetic energy of cold neutrons. Cold 
neutrons are therefore used in UCN production exper- 
iments as they can efficiently couple to the condensed 
matter excitation modes and exploit the available energy 
exchange channels. Our goal is to transform as many 
incident cold neutrons into UCN as possible. Neutron 
downscattering through excitations in a bulk material is 
known to be the most efficient mechanism by far for this 
process. 

In a solid target, the anisotropy of the lattice struc- 
ture manifests itself in an orientation-dependent disper- 
sion relation, as illustrated in FigEl As a result, for every 
neutron momentum transfer k, there exists a continuum 
of energy modes available for excitation. A large fraction 
of the incident cold neutron spectrum can therefore be 
used for UCN production. This is a potential advantage 
to the case of neutron scattering in super-fluid 4 He (the 
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FIG. 3: The density of states of magnon excitation modes in 
solid a-02, dark histogram for the acoustic mode and light 
for the optical mode. 



first proposed UCN source) , where only a narrow band of 
incident neutrons with energy centered at 11 K can pro- 
duce UCN. This results from the degenerate dispersion 
curves of phonon excitations in the super-fluid, which is 
an isotropic media. 

When sampled uniformly over the 3-D magnon mo- 
mentum space <f, the density of magnon states Z(E) can 
be histogramed over the energy scale, according to the 
allowed energy given by the dispersion relation Eq.@, 
i.e., 



J2 = JdEZ(E). 



(7) 



The density of magnon states in S-O2 is plotted in Fig. 13 
A distinct peak at 10 meV results from the flattening of 
the dispersion curve approaching the Brillouin zone (BZ) 
boundaries of n x . One would also expect a 20 meV peak 
corresponding to the BZ boundary of n y , however, this 
peak is less significant because of the smaller flattening 
region due to the higher periodicity of n y . Note that 
there are no magnon modes with energies below Eo = 
0.7 meV. The absence of low energy magnon modes is 
particularly important for UCN upscattering events, and 
we will return to this point when discussing UCN loss in 
Sec. III. 



C. Scattering Law in AFM crystals 

The reference by Lovesey^ contains detailed calcula- 
tions liking the scattering matrix elements to the scatter- 
ing law in Eq. (J2J) ■ The scattering law for neutorn inelastic 
scattering from a Heisenberg anti-ferromagnet with sin- 



gle magnon creation(+) and annihilation(— ) is given as This leads to 
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where F (re) is the electron spin form factor of an oxygen 
molecule, g = 2 is the electron gyro-magnetic ratio, k x 
is the unit neutron momentum vector projection on the 
magnetization direction, N is the total number of scat- 
terers, vq is the unit volume occupied per scatterer, and n 
is the magnon occupation number. For the zeroth order 
approximation^ we use the spherical Bessel function 15 
of order zero, jo(rei?/2), for the electron spin form fac- 
tor F(re). The internuclear distance R is 1.21 Aw& The 
additional terms u q and v q are due to the two sets of 
spin waves in the opposite sublattices, and they can be 
reduced to 



{ Uq + (-i) m v q y 
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with m = for a nuclear peak, and m = 1 for a magnetic 
peak. 

Summing over the allowed momentum transfer q, 
Eq. (JSJ can be further simplified to terms involving a dis- 
crete sum over allowed reciprocal lattice momentum r: 
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D. Evaluation of Neutron Coherent Scattering 
from Polycrystal 

We are now ready to perform the cross-section calcula- 
tions with the scattering law S(re, to) for S-O2. We present 
here our algorithm for calculating the pure coherent scat- 
tering cross-sections. We are interested only in the total 
cross-section which, according to Eq. J2J , is 



a tot {E u iU) =rl J dE f f tf%^S(re, 



co). (11) 



To consolidate the number of variables, we apply a co- 
ordinate transformation (CI, Ef) — > (re, to) — ((f), re, uj) 
with 
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Here, the neutron energy transfer u> is also the magnon 
energy gain, and re represents the magnitude of neutron 
momentum transfer. The total scattering cross-section 
becomes 
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a tot (E i ,n i ) = r 2 /// dudKtty-^'-Z-Sfcw), (13) 
which can be further reduced to 
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In a typical material with a periodic lattice structure, 
the dispersion relation contains the maximum energy at 
the Brillouin zone boundary, i.e., uj max = Ebz- In the 
case of single magnon exchange, this defines the lower 
bound (energy gain) of the energy transfer integral. The 
upper bound (energy loss) can be up to E BZ or the initial 
energy of the neutron. The lower, Ki ower , and upper, 
K upper , bounds in the momentum integral correspond to 
forward and backward scattering events, respectively: for 
a given magnitude of energy transfer and initial neutron 
energy, 
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The bounds 
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(y/Ei ~0J+ y/Hi) (17) 



are plotted, together with a contour plot of the scattering 
law, in Fig0] The area bounded by these two curves 
defines the integration area. 

Note that in the reduction Ea. l|13Jl to Ea. (|14|l . we have 
averaged the scattering law S(re, u) over the incident neu- 
tron angle £7,;, i.e., 



(18) 



<S(re, uj) = — / dfliS(H, u). 



The reduced scattering law, S(k,u>), is made isotropic, so 
that the double scalar variable integration in Eq. (|14|) be- 
comes straightforward. Note that even though it results 
in the same form for the final integral, the treatment 
presented here does not make use of the "incoherent ap- 
proximation" widely used in other cross-section calcula- 
tion packages^ The incoherent approximation relaxes 
the momentum conservation essential to coherent pro- 
cesses, and therefore it is subject to errors in the low 
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FIG. 4: A contour plot of the neutron scattering law S(k,lu) 
in a 5 K solid C1-O2: The upper half plane(a; > — ► E, > 
Ef) covers neutron downscatterings, and the u < covers 
neutron upscatterings. The area of integration is bounded by 
the dotted lines for incident neutrons with an initial energy 
Ei=0.03 eV, and solid lines, dashed lines and dash-dotted 
lines, for E;=0.01 eV, 3xl0" 3 eV and 10~ 5 eV, respectively. 



energy limit. Here, we arrive at Eq. l|14|l after integrat- 
ing the scattering law §(/?, to) over the incident angle, 
and this is suitable for the case of neutron scattering off 
polycrystals. 

Ea. i|14|) is evaluated numerically using a Monte-Carlo 
technique. We construct a table of values for the scatter- 
ing law S(k,lu) using the dispersion relation uj(H) given 
in Eq.©. With a given k, the vector k= (k sm. 9 cos <fi, 
k sin 9 sin <f>, k cos 9) is specified for a random choice of ori- 
entation angles (9,<fi); the dispersion relation u)(k) gives 
the corresponding energy of excitation, lu. Momentum 
transfers larger than the BZ momentum are mapped back 
to the first BZ, according to Ea. (|10fl . via unique Bloch 
translations. A numerical value of S(Ki,uij) is then cal- 
culated and added to a running total for each bin (i,j) in 
the parameter space. This process is repeated, stepping 
through the parameter space according to zAk = m and 
jAlu <u><(j + 1)Alu. Finally, we implement the angu- 
lar average by repeating the Monte Carlo random angle 
assignment for each k, and histogram the corresponding 
S(k,uj) values. 

Combining all the ideas above, we proceed to integrate 
Eq. (|14|l numerically by summing over all kinematically 
accessible points on a grid in the (u>, q) phase space. In 
the numerical approximation, the neutron downscatter- 



FIG. 5: Neutron scattering law S(k,u>) in a 5 K solid a-02'- 
S{k,u) is plotted as a function of energy transfer lu for mo- 
mentum transfer Q=0.1, 0.3, 0.5, 2.1, 4.1 and 6.1 A . Dis- 
tinct energy dependences are shown for lu > downscattering 
and lu < upscattering parts. 



ing cross-section is given by 
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For the results presented here, we constructed a table 
of S(k, lu) values on a grid of 500 x 500 points, equally 
spaced by An = 0.02A and Alj = 0.04 meV. The scat- 
tering law in a 5 K sample of solid oxygen is presented as 
a contour plot in Fig. 0] and the detailed dependence of 
S(k, u>) on energy and momentum transfer is illustrated 
in Fig. [S] Notice that in both figures, at low k, there 
exist two dominant contributions at 0.7 meV and 3 meV, 
which correspond to the self-energy terms involving D 
and D' of the acoustic and optical modes, respectively. 
As k increases, the peak in the scattering law shifts to a 
higher value centered at 10 meV, which is clearly the BZ 
boundary energy, as discussed before. 

Double integration of S(k,ui) using Ea. l|19|) and 
Eq. (l20|) then produce the downscattering and upscatter- 
ing cross sections for neutrons with a given initial energy. 
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FIG. 6: Cross sections as a function of the incident neutron 
energy: Solid curves: neutron upscattering from solid a-02 
at temperatures of 1 K, 2 K, 3 K and 5 K with annihilation 
of a single magnon. Circled curves: neutron downscattering 
via creation of one magnon from solid oxygen of 1 K, 2 K, 
3 K and 5 K. Note that the downscattering cross section is 
nearly independent of the solid temperature. Dashed curve: 
neutron-phonon scattering in solid oxygen. Dotted curve: nu- 
clear absorption cross-section per 16 02 molecule. 



The sum of the two gives the total cross sections. The 
resulting neutron-magnon scattering total cross sections 
in solid oxygen at various temperatures below Eo are pre- 
sented in Fig. 



E. Phonon Scattering in S-O2 

Neutron-phonon inelastic scattering is also present in- 
side solid a-02- The scattering law for single phonon 
creation(+) and annihilation (—) is^ 



§ ± (k,w) = 



1 



2M, 



,-VW{H) 



Z)W 



VQ 



E 

3=1,2,3 



(21) 



E '"^( ( f 2 (" + \ ± ^(^.- T M<* 3 (? T H r). 

q.T J \"J 

The neutron-phonon scattering cross-sections are cal- 
culated using LEAPR, one of the subroutines in the 
NJOY package developed in the Los Alamos National 
Laboratory^ The incoherent approximation is used for 
the oxygen nucleus O 16 , even though it is a pure coherent 
scatterer. While we still expect a discrepancy at low en- 
ergies, we expect the incoherent approximation to yield 
results differing from the actual phonon cross section by 
at most a small numerical factor—. An iterative expan- 
sion is used to include creation and annihilation of up 
to 100 phonons. The cross section for neutron-phonon 



scattering is shown along with those for neutron-magnon 
scattering in FigEJ 

Note that the above formula is similar to Eq. JHJ , with 
the important distinction of an additional inverse depen- 
dence on the scatterer mass M s . This indicates that 
phonon production is suppressed due to the smaller am- 
plitudes associated with heavier nuclear recoil. Con- 
versely, magnon excitation takes place via neutron- 
electron interactions, and thus does not have the mass 
suppression. To be more precise, we will compare the 
prefactors and the energy terms in Eq.JHJ and Ea. f21|) 
with the corresponding scattering length: for the magnon 
part, one finds 
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for the phonon part, one has to pay attention to the 
detail that there are in fact two nuclei in one molecular 
scatterer, which contribute to the scattering amplitude 
coherently. Therefore, the relevant term is 
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where A is the atomic number. As the intrinsic nu- 
clear coherent scattering length (b c — 5.78 fm for the 
16 O nucleus )il is comparable to the magnetic scattering 
length (ro = 5.4 fm), the ratio of the neutron-phonon 
cross section to the neutron-magnon scattering cross sec- 
tion would naively be on the order of 1/A = 1/16, assum- 
ing similar kinematics. The results of numerical integra- 
tions are presented in FigH3 which show that the neutron- 
phonon scattering is overall weaker than the neutron- 
magnon scattering by about one order of magnitude for 
incident neutrons with energies less than 10 meV. This 
is consistent with the above analysis. 



III. COMPARISONS OF PERFORMANCE 
BETWEEN S-0 2 AND S-D 2 

A. UCN production 

For the purpose of constructing a practical UCN 
source, we are interested not in the total cross section, 
but in the scattering events that cover the UCN phase 
space. To estimate the UCN production rate given a flux 
of incident neutrons with a broad spectrum, we carry out 
an integration of the differential cross section weighted by 
the cold neutron spectrum, i.e., 
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gen, by contrast, the incoherent elastic scattering cross- 
section vanishes, so that the elastic mean free path is in 
principle infinite. Thus the fiducial volume of the solid 
oxygen source is expected to be limited solely by the 
neutron absorption length of 380 cm. Therefore, an im- 
provement in source volume by a factor of up to 



380 



50 fl 



d=l, 2, 3 



(24) 



over solid deuterium is theoretically possible, though this 
will be limited by the maximum practical size of a solid 
cryogenic source. Here, d is the number of dimensions in 
which the source volume can expand. The tendency of 
micro-crystals to form in S-C^r^^ however, may limit 
the effectiveness of a large UCN source. 



Cold Neutron Temperature (K) 



FIG. 7: UCN density vs. energy of the incident cold neutron 
flux characterized by the temperature a Maxwell-Boltzmann 
spectrum. The UCN production is optimum with a flux of 
12 K cold neutrons in S-O2, whereas it is optimum with a 
flux of 30 K cold neutrons in S-D2. 



where n is the density of the scatterers. We assume a 
Maxwell-Boltzmann energy distribution for the incident 
neutrons. The results are plotted in FigEI as a func- 
tion of the incident cold neutron temperature. With a 
typical flux of 30 K cold neutrons on solid oxygen, we 
predict a UCN production rate of 2.2 x 10 _8 $o, where 
$0 is the total incident cold neutron flux. As shown 
in the plot, the optimum cold neutron temperature for 
UCN production in solid oxygen is about 12 K, which 
is roughly half the optimum cold neutron temperature 
for a S-D2 source. This gives a UCN production rate 
of P UC n — 3.5 x 10~ 8 $0i a 60% improvement over the 
30 K cold neutrons. For comparisons, the UCN produc- 
tion rates in S-D2 are 1.95 x 10 _8 $o (in pure ortho S-D2) 
and 1.74 x 10~ 8 $ (in pure para S-D 2 ) with 30 K cold 
neutrons. 

The maximum achievable UCN density is the produc- 
tion rate times the lifetime r, i.e., 



= Pu 



X T. 



(23) 



For S-O2 at 2 K, if there are no loss mechanisms other 
than nuclear absorption and magnon upscattering, r = 
375 ms. By comparison, the longest lifetimes achievable 
in S-D2 are approximately 40 ms^ with the relative re- 
duction due primarily to the larger nuclear absorption 
cross-section, hydrogen impurities and finite concentra- 
tion of para-molecules. As a result, an improvement in 
UCN density of about an order of magnitude is expected 
in solid oxygen relative to solid deuterium. 

Finally, the maximum achievable UCN flux is also de- 
termined by the fiducial source volume. For a solid deu- 
terium source, the fiducial volume is limited by the in- 
coherent elastic scattering length of 8 cm4 In solid oxy- 



B. UCN loss 

1. Thermal Upscattering 

The 2-D symmetry of the AFM magnons in S-O2 re- 
sults in a fundamentally different dispersion curve from 
that of phonons. As it turns out, it takes a finite mini- 
mum energy E of about 8 K to excite a magnon. For 
a solid oxygen sample of volume V with a temperature 
below 8K, the total number of magnons is given by 



N = J2n t S 



V 



(2ny j Eo 

e -8K/T 



Z(e)de 



1 



e/k B T _ 1 



(25) 



That is, the number is exponentially suppressed as the 
solid coolsi2i For temperatures lower than 8 K, there are 
no low energy magnon levels available to be thermally 
populated, and only a small population of magnons with 
energy above Eo exists. This results in a partial freeze- 
out of the UCN upscattering mechanism. The resulting 
exponential factor, e~ 8K ' T , of the magnon annihilation 
process results in a stronger temperature quenching ef- 
fect on UCN upscattering than that of the typical phonon 
absorption mechanism. Note that the loss rate due to 
magnon upscattering is reduced to levels comparable to 
the nuclear absorption rate when the solid is cooled down 
to 2.1 K (see Fig. EJ. This sets the operational temper- 
ature of a solid oxygen UCN source, as further cooling 
should increase the lifetime by no more than a factor of 
two. 



2. Para-molecule Upscattering 

Natural oxygen is over 99% 16 0, which has a nuclear 
spin of zero. The molecular rotational spectrum has only 
one set of allowed states, and therefore 16 02 molecules 
can not be left in long-lived rotational states at low tem- 
perature. The interaction of UCN with the rotational 



states of D2, however, is a major source of loss in S-D2 
sources, and must be controlled by minimizing the para- 
D2 concentrationi^fliSi UCN in S-O2 are exempt from 
such a problem. 



tering. Oxygen appears to be a promising candidate, 
with strengths and weaknesses which complement those 
of solid deuterium. 



3. Losses specific to S-O2 

While nuclear absorption and magnon upscattering 
losses are theoretically well understood, this is not the 
case for the effects of polycrystal formation and various 
lattice defects. Previous experience with solid deuterium 
indicates that finite crystal effects do not introduce ad- 
ditional measurable scatterings, however, this statement 
might not apply to solid oxygen because of its very dif- 
ferent thermal properties. This can only be answered by 
experiments. 

In high radiation fields typical for neutron sources, 
ozone formation^ in oxygen solid may be a source of 
UCN loss, as the increase in mass fluctuations produces 
additional incoherent scatterings^ We also expect more 
7-ray heating in oxygen because of its higher Z relative 
to solid deuterium. This, together with a lower thermal 
conductivity^ may increase the demand on cryogenic en- 
gineering to maintain the solid oxygen at 2 K. We have 
also estimated the effect of additional hyperfine transi- 
tion energy in the S-O2 crystal^ on UCN loss. While we 
expect no significant rate of UCN upscattering from this 
mechanism, definitive measurements of the UCN lifetime 
in solid oxygen are required to place quantitative limits 
on this and other loss mechanisms. 



IV. CONCLUSIONS 

We have reviewed the potential of solid oxygen as a 
UCN source, by analyzing inelastic neutron-magnon scat- 



The physics of UCN production in solid oxygen, in- 
volving magnon (spin wave) exchanges, is fundamentally 
different from the well-known phonon mechanism in solid 
deuterium. In solid oxygen, the neutrons are down- 
scattered through coupling to the magnons in the anti- 
ferromagnetic a-phase of the solid, the strength of which 
is comparable to the nuclear scattering (phonon) in solid 
deuterium. This mechanism has several advantages. For 
one, the absence of magnon states at low energy implies 
a greatly reduced UCN upscattering probability, relative 
to solid deuterium. Second, the smaller nuclear absorp- 
tion cross-section of oxygen leads to a much longer UCN 
lifetime in this material. Also, the UCN loss due to the 
absorption of the excess energy from para-molecules (a 
major source of loss in solid deuterium) is absent in oxy- 
gen. 



Preliminary calculations for solid oxygen indicate that 
the UCN density achievable should be about an order 
of magnitude greater than for solid deuterium. A fac- 
tor of 10 gain in source volume should be readily feasible 
by taking advantage of the smaller neutron absorption 
probability and the infinite elastic scattering mean free 
path of UCN in solid oxygen. Together, these proper- 
ties of solid oxygen suggest that a hundredfold increase 
in the UCN flux relative to what is available from solid 
deuterium should be possible, given the same incident 
neutron flux. The optimal operational temperature is 2 
K, somewhat lower than that of solid deuterium (5 K). 
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